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Abstract 

We discuss the renormalization of the Newton constant due to fields of various spin s. 
We first briefly review the cases of s = 0, 1/2, 1, 3/2 already discussed in the literature 
and notice the appearance of the well-known contact terms for the vector bosons. We 
then extend this discussion of the contact terms to massive vector fields, p-forms and 
to the case of spin s = 2 particles (gravitons). We observe that, in general, the contact 
terms originate from the fields which mediate the interactions (such as vector gauge bosons 
and gravitons). We then discuss entanglement entropy and the conical entropy and their 
relation to the renormalized Newton constant. We address the puzzle of the non-analytic 
terms due to fields of spin s = 2 and suggest that the resolution of this puzzle comes from 
the non-equivalence of the orbifold and n-fold cover constructions which are used in the 
entropy calculations. Finally, we propose a mechanism by which the Bekenstein-Hawking 
entropy is identified with entanglement entropy in any theory which includes both matter 
fields and the mediators of interactions (vector gauge bosons and gravitons). 
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1 Introduction 


The Newton constant not only defines the strength of gravitational interaction in a low energy 
regime but also determines the value of entropy associated with the black hole horizons [1]. Both 
aspects arise in a semi-classical approximation in which the graviton loops [2] do not play an 
important role and, as it seems, can be ignored. The statistical aspect of horizons, however, is 
believed to be a window through which one may try to see the fundamental degrees of freedom of 
Quantum Gravity otherwise unaccessible. In the recent years there has been proposed a number 
of approaches that explore this window in various directions 0. a. 0. On the other hand, 
it seems plausible that, in a quantum theory, the observable Newton constant and the entropy 
should arise simultaneously within the same mechanism which deals with the observable low 
energy particles, irrespectively to the details of the hypothetical more fundamental description. 

The resolution of the problem in this formulation boils down to the analysis of the renormal¬ 
ization of the Newton constant and the appropriate calculation of the entropy [6]. The entropy 
most suitable to this purpose is the entanglement entropy [7j defined for a given surface as en¬ 
tropy of ignorance of the modes residing inside the surface, for the earlier works in this direction 
see [8] and more recently 0, d, mi, n. This definition is very much suitable to describe the 
situation with the principle impossibility to gain information on that physical modes which prop¬ 
agate inside a black hole (for a review see |T3], [H]). The departure point in this direction is the 
observation that the entanglement entropy, if combined with the classical Bekenstein-Hawking 
entropy, can be expressed in terms of the renormalized Newton constant 0, IE], 03, S3, |I3|. 
If it worked for all fields this would allow us to relate directly the two aspects, the generation of 
the observable Newton constant and the statistically meaningful Bekenstein-Hawking entropy. 
It was observed, however, very soon that this does not work for all fields. Namely, it does not 
work for a scalar field non-minimally coupled to gravity [19] and for the gauge vector fields [21]. 
Kabat has observed that the difference is due to the presence of the so-called contact terms 
which do not appear to have a statistical meaning. The issue of these terms has been actively 
discussed in the literature [23]. These terms are often considered as obstruction to the success¬ 
ful realization of the program of identifying the Bekenstein-Hawking entropy as entanglement 
entropy. 

In the present paper we extend the discussion of the contact terms to other fields: massive 
vector fields, p-forms and particles of spin s = 2 (gravitons). The result for the contact terms 
due to gravitons was already given (without details) in [13]. Here we give a detail analysis 
of the linear perturbations in the Einstein-Hilbert action in arbitrary dimension d > 4. We 
decompose the gravitational perturbations of the Einstein-Hilbert action on the trace (conformal 
mode) and the traceless part and show that in generic dimension d > 4 the conformal mode 
becomes non-minimally coupled so that the coupling vanishes only if d — 4. Moreover, we show 
that the relevant operator which governs the traceless perturbations is a curvature-dependent 
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modification of the Lichnerowich operator. The complete treatment of all perturbations allows 
us to determine the contact terms due to gravitons. In general, we notice that the presence of 
the contact terms is due to the fields which mediate the interactions and thus are due to the 
presence of the (exact or broken) gauge symmetry. 

We then suggest that, considering all fields that appear in a low energy regime, there exists 
a particular choice of the otherwise arbitrary bare gravitational constant that the Bekenstcin- 
Hawking entropy, defined in terms of the renormalized Newton constant, is identical to the 
entanglement entropy. The UV cut-off needed for the proper definition of the entanglement 
across the horizon then is automatically set by the Planckian scale. In this mechanism the 
two pieces in the Bekenstein-Hawking entropy which do not have a statistical meaning, the one 
which is defined by the bare Newton constant and the other due to the presence of the contact 
terms in the renormalized Newton constant, mutually cancel so that what is left is genuine 
entanglement entropy. 


2 Newton constant in quantum action 


We start with reminding a general set up for the calculation of the Newton constant which 
appears in the quantum effective action. The effective action of a field of spin s can be written 
as 


W { .) = 


(-) 


2s 


r ^Tre- A(3) 

,2 r 


( 2 . 1 ) 


The second order covariant operators acting on the spin- s held can be represented in the fol¬ 
lowing general form 


A (s) 


—V 2 <5 


AB 


I y( s ) 
+ A AB 


( 2 . 2 ) 


where the matrix X( s ) depends on the chosen representation of the quantum held and it is linear 
in the Riemann tensor. The UV divergences of the action (12.11) are analyzed using the small 
proper time r for the heat kernel of this operator 


Tre -sA(s) (r) 



y^ y a n T 

n =0 


(2.3) 


The induced Newton constant is determined by the coefficient ai in this expansion. It has the 
general form 

ai s) = f d d x^(^-R - TUf«), (2.4) 

where D s (d) = Saa is the dimension of the representation of spin s. For s < 2 spins we have 


D s=0 (d) = 1, D s=1/2 (d) = 2^ , D s=1 (d) = d , 

D s=3/2 (d) = d2W* , D s=2 {d) = (rf ~ 1} 2 (rf + 2) . (2.5) 
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D s (d ) can be interpreted as the number of off-shell degrees of freedom, for s = 2 it counts 
number of degrees of freedom of a symmetric traceless matrix. 

The effective actions is expanded in powers of the curvature. The first term in this expansion 


is 


K! - 


16n G(e) 


d d Xy/gR = 


(- 1 ) 


2s 


» 


(An) 2 (d — 2) 


cd- 2 


( 2 . 6 ) 


where G(e) is the induced Newton constant, e is the UV cut-off in the theory. 

In what follows we shall give the calculation of the induced Newton constant for fields of 
various spin s. The result can be represented in the following form 


1 


1 


1 


cd- 2 


V s (d) 

6 


-C(*)(d) 


(2.7) 


4G(e) (4n-) d 2 2 (d — 2) 

where T> s (d) is the number of on-shell degrees of freedom, each fermionic degree of freedom is 
counted with the weight 1/2, 


2 [d/2] 

T> s=0 (d) = 1, T> s= i/2(d) = ——— , T> s= i(d) = d — 2 

oW 2 } d(d — 3l 

V p (d) = C p d _ 2 , V s=3/2 (d ) = (d - 3) — , V, =2 (d) = . (2.8) 

where p stands for p-forms. The quantity c^(d) which appears in (I2.7jl is what we shall call 
the contact term for field of spin s. This quantity is non-zero only for the bosonic gauge fields: 
s = 1 vector bosons, s = 2 gravitons, p-forms, 

c { ’Hd) = ?, = 1, c<?>(d) = C£l , c‘- 2| (d) = r LlA±± , (2.9) 

where ^ is non-minimal coupling of scalar field to gravity. The cases s = 0, 1/2, 1 were already 
considered in the literature. We shall give below the detailed analysis of the case of s = 2 spin 
and of the p—forms. 

We start with some well-known examples (see eu, m, [Z 2 ], my 


3 Fields of spin s < 2 

3.1 Non-minimally coupled scalar field, s = 0 

A non-minimally coupled scalar field is described by the operator 

A( s= °) = - V 2 + £R , (3.1) 

where £ is a parameter, £ = £(<f) = -yyvjy corresponds to the conformal case. The contact term 
in this case is 

c (s=o) (d) = £ ■ (3.2) 
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3.2 Dirac fermion field, s = 1/2 

The partition function of the Dirac fermion, Z\j 2 = det 1 ^ 2 A^ s=1//2 ^ is expressed in terms of the 
squared Dirac operator. The latter takes the form (12.21) with matrix 

At ' /2) = \rSab , (3.3) 

where R is the Ricci scalar and (A, B) are the spinor indexes. Thus we have that 

(-l) 2 *a<- 1/2 > = -D„ 1/2 (d)( i )Jr= P - 1 6 /2(d) J R . (3.4) 

The comparison with (I2.7[) shows that the contact term vanishes in this case 

c (s =i/ 2 )(d) =0. (3.5) 


3.3 Gravitino, s = 3/2 

The partition function, including the Faddeev-Popov and the Kallosh-Nelson ghosts contribu¬ 
tions, is Z 3/2 = det 1/2 A( 3 / 2 ) def 3/2 A^ 2 ) , see for instance 20]. The operator A^ 3 / 2 ) takes the 
form (1272|) with 

X A&hv = -^R^AB9fiu - -jR^ivcxP^l^AB , 

TrX^ = ^D 3/2 (d) J R. (3.6) 

The heat kernel coefficient of this operator reads 

A = D 3/2 (d)(\ ~ j) J R = ~D vl (d) jR (3.7) 

so that, taking into account (13. 4 j) . we find that the appropriate heat kernel coefficient in the 
partition function is 

- ai = -a 3 / 2 + 3 a 1 / 2 = -V 3/2 (d) j R, (3.8) 

where V s = 3 / 2 (d ) = (d — 3) i s the number of on-shell degrees of freedom for spin s = 3/2. 
Notice that it vanishes in d = 3 as it should be by supersymmetry! 3 ]. The comparison with (12.71) 
shows that the contact term vanishes in this case too. 

1 We thank M. Duff for this remark. 
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3.4 Gauge field, s = 1 

The case of gauge vector bosons has been actively considered in the literature. Here we briefly 
summarize the analysis for the Newton constant. The partition function of a gauge boson is 

Z (1) = det _1/2 A^ ) det A (0) , (3.9) 

where Ad) is the operator of the form (12.2p with 

XjS = +R IM ,, Tr xV=R. (3.10) 

The heat kernel coefficient of this operator is 

4 1 ’= (g-i) J r- 

The ghost operator in (13. 9 j) is the minimal scalar (s = 0) operator 
the ghost contribution one has that 

a x = - 2ao S-0) = (" 1 - jj) J R= (^s=i( rf ) “ 

The contact term thus is 

c (a=1) (d) = l. (3.13) 

This is the famous Rabat’s contact term [2TJ (see also [22])- We see that the induced Newton 
constant due to the spin s = 1 field is not a priori positive. Indeed, it is negative if d < 8, 
vanishing if d = 8 and positive if d > 8. 


(3.11) 

. Thus, taking into account 
-1 )Jr- (3.12) 


3.5 Massive fields 


In the case of scalars or Dirac fermions it is rather straightforward to add mass into consideration. 
One simply adds a mass term to the second order operator A Sim = A s + m 2 , s = 0, 1/2. 
Respectively, one modifies the heat kernel, 

Tr e ~ TAs ’ m = Tre~ rAs e ~ Tm2 , (3.14) 


where Tre tAs has the decomposition (12.31) . The induced Newton constant then is given by the 
integral 


1 

4 G{e) 


V £ 
12 



dr 

(4 Trr) d / 2 


e - m V 


s = 0, 1/2 


(3.15) 


To leading order in e this gives us the same behavior as before for the massless fields. On the 
other hand, there appear some sub-leading terms as function of m 2 e 2 . 


7 







The analysis is more subtle for a massive vector field. The reason is that the gauge invariance 
is broken by the presence of the mass term that effectively modifies the number of the physical 
degrees of freedom. The analysis made in [26J shows that the partition function of a massive 
vector field is 

Z s= i )Tra = det -1//2 (A^ + m 2 ) det 1,/2 (A® + m 2 ). (3.16) 


The respective contribution to the Newton constant is 


1 

4 W) 


1 ( 'E > s=l,m _ (s=l,m)\ 

2 \ 6 ) 


r dT c - m2r 

J e 2 (47rr) d / 2 


where 


(3.17) 


V s=1 . m (d) = d- 1, c<- =1 > m >(d) = l. 


(3.18) 


We see that the structure of the induced Newton constant remains the same as in the massless 
case. The physical degrees of freedom (their number is modified if compared to the massless 
case) contribute positively with the weight 1/6 while the contact term contributes negatively to 
the gravitational constant. Notice, that the contact term in (j3.18[) is the same as in the massless 
case. 

Generalizing (13. 1 7ft to fields of various spin s and mass m s we have 


1 


'Em - r dr c 

^ 2 \ 6 J Je 2 ( 4:TlT ) d/l2 


for the complete contribution to the Newton constant. 


(3.19) 


4 Gauge £>-forms 

The other example of lower spin fields with gauge symmetry are the gauge p-forms, 

uUi, i dx 11 A .. A dx lp , 

p\ 

where is a totally antisymmetric tensor of rank p. The action 

W (p) = ^J *du {p) A duj {p) 

is invariant under the gauge transformations, 

Jp) u (p) + da (P~ i) ) 


(4.1) 


(4.2) 


(4.3) 


where a^ p ^ is a (p — 1)-differential form. The p-forrn is thus a generalization of the Maxwell 
gauge held considered in the previous section. The fixing of the gauge symmetry (I4.3P requires 














introduction of a set of the Faddeev-Popov ghosts, the differential forms of rank p—k , k = 1,2,.. 
some of them considered as bosonic fields and the others as fields of fermionic statistics. The 
resultant partition function then takes the form (see for instance [28] ) 


^) = n ( detA p _ fc )(- 1 ) fe+1 ^ i , (4.4) 

k =0 

where A fc = d k ~ 1 6 k + S k+1 d k is the Beltrami-Laplace operator, a generalization of the Laplacian 
for the k -forms, which is expressed in local coordinates by the Weitzenbock formula [29] 

k 


(A^ ..ifc 




ni\ ..ii/ ..ik 


+ 2^(-l Y +v R r 




(4.5) 


v=\ 


[1<V 


It takes the form —V 2 + X^ k \ where the exact form of X^ can be deduced from (14.51) . The 
calculation shows that TrAiW = R, TrA^ 2 ) = (d — 2 )R, TrX ® = ( d ~ 2 K d ~ 3 ) r The natural 
extrapolation to arbitrary k is TrA^ fe ) = = C^zl R- In dimension d = k the k- 

form is dual to 0-form and thus the operator A*, becomes the usual Laplace operator acting on 
scalar functions. Therefore Tr = 0 if d = k for k > 1. 

The heat kernel coefficient corresponding to the product of determinants in (14.41) is a sum 
of coefficients for operators (14.5ft of various rank, 

«?’ ) = E(*:+l)(-l)'‘« 1 (Ap- l ). (4.6) 

k =0 

For the heat kernel coefficient of operator A*, we have the formula (12.41) 

ar(A*) = f d d xj-g - Tr , D k (d) = C k , TrX™ = C^R, (4.7) 

where we take into account that the number of independent components of an antisymmetric 
tensor of rank k in d spacetime dimensions is C\. 

Now using the identity 


XTr‘(-i)‘(*: + i) = cy 2 


the sum in 


k =o 


can be evaluated and we fold 


thatj 2 ! 


o?> = [ d d x ys ( 


V 6 


-c ip \d)^j R, 


where T> p (d) = C p _ 2 is number of on-shell components of 71 -form and 


(4.8) 


(4.9) 


c M (d) = Ctl 

2 As we have learned this result was earlier obtained in [3D]. 


(4.10) 
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is the corresponding contact term. This expression generalizes Kabat’s contact term (j3. 13j) to 
arbitrary gauge p-form. Notice that T> p (d) as function of dimension d is a polynome of degree 
p while c {p) (d) is polynome of degree (p — 1 ). 

Using the relation p(d — 2 — p)C%_ 2 = (d — 2)(d — Z)C^Z\ we can re-write (j4.9|) as 

a?' = CZ / -m. (4.11) 

This shows that the corresponding contribution to the induced Newton constant is mostly 
negative for p < d — 3. The only case (for p > 1) when it exactly vanishes is d = 8, p — 5. 

5 Graviton, s = 2 

Our next example is the graviton. This, obviously very important case, is less studied in the 
literature. Therefore, we shall give as much details as possible. The first important question 
here is what the appropriate differential operator which describes the quantum dynamics of 
the graviton? In the literature one typically uses the Lichnerowicz operator. However, the 
Lichnerowicz operator is not the one who governs the dynamics of the quadratic fluctuations in 
the Einstein action. In order to find the form of the appropriate operator we shall study the 
quadratic term in the Einstein action 

WeH = _ 167tG / ddx ^ R l 5 - 1 ) 

with respect to perturbation of the metric, g = gffi + \/l 67 r ~Gh tlv , where gjS is the background 
metric and h pi/ is the perturbation. We stress that we want to see what Newton constant is 
induced by these fluctuations. Therefore, the background metric is not supposed to be con¬ 
strained by any field equations. In what follows we shall skip the script (0) for the background 
metric. 

5.1 Linear perturbations and the Lichnerowicz operator 

To start with let us consider the linear perturbation of the Ricci tensor, 

SiRpv = - V 2 /r Mi , - V^VjTr) 

= - \s‘h) + - \s°h) + 1a l / v , (5.2) 

where in the second line we commuted the covariant derivatives and A l is the Lichnerowicz 
operator 

A L h Q/ 5 = —V 2 h a p + Rfhtf + Rq h m - 2R' 1 Q l ' ff h lw . (5.3) 
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(5.4) 


Assume that the linear perturbation satisfies the gauge condition 

V“(h a/3 - ^g a ph) = 0 , h = . 

one finds that the linear perturbation of the Ricci tensor 

Rfiu 2 ^Lhfj,u . (5.5) 

is governed by the Lichnerowicz operator. In fact, the Ricci tensor is the only tensor the 
small deformation of which is determined by the Lichnerowicz operator without any curvature 
dependent modifications. 

Respectively, the linear perturbation of the Einstein tensor G^ v = R^ w — | Rg ^ reads 

Si GV = - ls?h) + - \s°h) - iv”V«(V - \g«f,h) 

+^A L h^ + \ 9lxv S7 2 h - ]ri^R + ^g^ha P R ap (5.6) 

In the gauge (15.4p the first line in (j5.6[) vanishes and we have 

Gfj U = ^A L h^ + ^g^ 2 h - hi^R + ^g P vh ap R rTp . (5.7) 

It is useful to introduce a traceless part of the perturbation, 0 M1/ = — ^g^h. Decomposing 

the Einstein tensor on the trace and traceless parts, 

8\G = g a ^8iG a p , SiG^p = 8\G a p — —g a p8±G ( 5 - 8 ) 

we find 

SiG = lAAv 2 /! + , 

$i G lp = ^ - ^({) a pR . (5.9) 

If the background spacetime is Einstein, i.e. R^ IU = Ag [W , then (f>p V R= 0 and the two modes 
(the trace and the traceless one) decouple. But even in this case the relevant operator, which 
governs the linear transverse traceless perturbations of the Einstein metric, in the Einstein equa¬ 
tions is operator (A l — R) which is a deformation of the Lichnerowicz operator by a curvature 
dependent term. We shall see in the next section that the deviation from A^ is even more 
important if we analyze the operator for the quadratic perturbations in the Einstein-Hilbert 
action. This consideration shows that the Lichnerowicz operator without any curvature depen¬ 
dent modifications appears only in the perturbation of the Ricci tensor. For other tensors, such 
as the Einstein tensor, the modifications due to some curvature dependent terms is necessary. 
Since the curvature dependent terms in A"^ s=2 ^ are essential for the structure of the induced 
Newton constant we conclude that it is not reasonable to consider the operator A l alone as the 
appropriate operator for the spin s = 2 particles. The curvature dependent modifications of 
A l should be considered that necessarily leads to modification in the induced Newton constant 
caused by the spin s = 2 particles. 
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5.2 Quadratic perturbations in Einstein-Hilbert action 

Let us now analyze what operator governs the quadratic perturbations in the Einstein-Hilbert 
action. For this we first consider the term quadratic in perturbations in the gravitational action 

(EZED, 

S 2 W EH = -J d d x^S7,{h^ - h)V a (h “ - l -5«h) 

- J d d x^h a0 A L h a p + hv 2 h - 2 h”G v a K - h 2 R + hh a pR a0 - hi a ph a0 Ri§.lO) 

where G a p = R a p — \g a pR is the Einstein tensor and we raise indexes with the help of the 
background metric g a0 . The term in the first line in (I5.10P originates from the first line in (15. 6p . 
As ’t Hooft and Veltman [2] have instructed us long ago it is convenient to add the following 
gauge fixing term 

W = 5 2 W eh + W„ , W„ = J d! J x^V IJ {h 111 ' - ^h)V a (hZ - l -5«h) (5.11) 

that cancels the term in the first line in (I5.10p . As in the previous section, it is useful to 
introduce the traceless perturbation = h a p — \gaph. Then the quadratic action (15.111) is 
expressed in terms of 0 a/ 3 and h 

w '=-[ - \kr)k 

In dimension d = 4 one reproduces the earlier results [31], [32]. 

Some remarks are in order. 

1. First of all it is important to note that the conformal mode h becomes non-minimally coupled 
if d 7 ^ 4. The corresponding operator takes the form (12.21) 

(-V 2 + XW ), X( h) = ——j——R. (5.13) 

LL 

This non-minimal coupling disappears if d — 4 but otherwise is generically present. To the best 
of our knowledge this fact is not known in the literature. 

2. There generically presents an off-diagonal coupling between and h. This coupling 
disappears only if d = 4. 

3. I 11 the traceless transverse sector we see the appearance of an extra term (j)"(R u a — |5“i?)0 a M . 
This extra term disappears in d — 4 if the background metric is Einstein so that the differential 
operator acting on the transverse traceless mode (j) a p is the Lichnerowicz operator. In dimensions 
d > 4 this terms does not disappear even if the metric is Einstein (provided the scalar curvature 
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is non-vanishing). Thus, in higher dimensions an important curvature dependent modification 
of the Lichnerowicz operator appears. 

The above observations show that d — 4 case is in many respects special. However, even in 
d — 4 for a generic (non-Einstein) background metric the operator in the transverse traceless 
sector is not the genuine Lichnerowicz operator. 


5.3 Newton constant due to graviton in d dimensions 

The quadratic action (15. lip should be supplemented with the ghost term which corresponds to 
our choice of the gauge fixing, 

W = w'+ W ghost , W ghost = i J d d x^-g - R^)xv ■ (5.14) 

The corresponding ghost contribution to ai is 

a f =^R- TrX gh , TvX gh = -R . (5.15) 


As it is seen from the structure of the action (15.121) the trace and traceless sectors should 
be considered together as a composite field (0^, h). The corresponding kinetic terrnj 3 ] however 
diagonalize in each sub-sector of 0and h. So that only the X -matrix has the off-diagonal 
(0/xi >h) components. These components however are not important when we compute the trace 
of the complete matrix X in the (0, h ) sector, 

Tr X^' h) = TvX tt + TrX (/l) , (5.16) 


where Tr A ^ is found in (j5. 13[) and matrix X T1 is dehned in the traceless sector. In this sector 
we have 


1 

2 


/ 


d d xXg (-yfv 2 + x™) <tr , 


(5.17) 


where 

C = \( S X + S K) - \s^S af ( 5 . 18 ) 

is the projector on the traceless part of a symmetric tensor, its trace = J "+ d " 2 = D s=2 (d) is 
the number of components of a symmetric traceless tensor of rank 2. Tensor X^ is defined 
as 


x: 


a/3,iiu X a p^ v ^gagX a p^ u g ^ d^^ LV ^ a ^ ,a pg ^2 CT p 5e/t <7 ^9 i (5.19) 


3 The kinetic term for the trace mode has a wrong sign. We use the standard prescription h —¥ ih. 
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where X a p^ u takes the form 

1 R 

Xap 

,fiu Rfiavp Rvafip 2 9Pli T ^a^i9/3i' “1“ R/3v9a[i T R/3fi9ceu) T ^ ^9jj.a9v6 T 9/j./39i'bfy 
The trace of (15. 19H is evaluated as 

-y-TT,a/3 _ Y a @ _ 1 Y n^ u n a 9 (A 911 

X a/3 X a p ^fxv,a09 9 ■ fO.Zij 


We find that 


X 


a/3 

afU 


= R 


d(d — 1 ) 


X^pg^g*? = R(d 


4) 


(5.22) 


and hence 


Tr.Y TT = X TT -% = R ( - * 2 + i). (5,23) 

The coefficient aq in the sector (</>, h) then 

a (M) = D s= 2 {d) + 1 ^ _ TrX ^,h), TrX (j>,h) = d(cl - 1) ^ 

6 2 

The above analysis in the traceless sector can be compared with the analysis for the Lich- 
nerowicz operator ( 15 . 31 ) . The matrix X in this case is 

x L ’Z = + \(K% + (/. ++ V, a « ft)) . ( 5 . 25 ) 

Note, that it is already traceless, X L,lx ”pg a P — 0 and X L '^g^ = 0. Then we find that 

TyX l = (d + 2 )R . ( 5 . 26 ) 


Curiously, (15.26P and (15. 23() . being generically quite different functions of dimension d, coincide 
if d — 4. This is due to the fact the difference of the matrices (15. 19j) and (15. 25(1 although 
non-vanishing has zero trace if d — 4. 

Returning to our analysis of the renormalization of the Newton constant due to the gravitorQ, 
we find that the complete coefficient a[ s- ^ is the sum of contributions of the ( 0 , h) sector and 


4 The earlier works, where a relevant analysis in dimensions d = 4 and higher was present, include [33] . [34] . 
[301 . In particular, in section 4.1 of |30 it is given a calculation similar to the one presented here. However, 
instead of separating the gravitational perturbations on the traceless and the trace parts as we do the authors 
of [30J deal with the operator acting on the complete gravitational tensor modes. The relevant operator is the 
one which appears in (15.1011 . (15.111) and it has the form (—PV 2 + X) where P and X are given in (4.5) of [501 . 
This operator factorizes on a product of P and (—V 2 I + P~ 1 X) so that only the latter produces an essential 
contribution to the determinant. The heat kernel coefficient is ai = Tr(P/67 — P -1 A'). One then finds that 
TV7 = d(d + l)/2 and TrP -1 X = Rd(d — l)/2 that exactly corresponds to our equation (15.241) . So that this 
way we again reproduce our results (15.241) and, taking into account the ghost contribution, (15.271) . The result 
given in (4.6) of [30] is, however, different from (15.241) by (d— 4)/2 and, as we believe, is incorrect. We thank 
A. Tseytlin for bringing to us his earlier relevant work and for the discussion of the relation between our result 
and the one given in |3U]. 
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of the ghost (15.150 . 


o (,.2) = a (4,,h) _ 2a ,k = 1 (Bj _ 2 + ! _ 2i)R _ f * 2 + * 

(d 2 -d + 4) 

6 2 


(^-4) 

d 


+ 2)R 


(5.27) 


where T> s = 2 (d ) = is number of on-shell degrees of freedom of graviton in dimension d. 

Notice, that in this expression the non-polynomial term 4/d of (I5.23f) cancels with —4/d coming 
from (15.131) . The contribution to the induced Newton constant thus takes the form (12.71) with 
the contact term 

/ i\ d ^ — d -\- 4 . . 

C(,=2)W = -2-' < 5 ' 28) 

This result (without details) is given in [f3j. Notice that the complete contribution to the 
Newton constant due to graviton is negative in any dimension d > 4. 


6 Some comments 

Let us make some comments. 

1. Our first observation is that the induced Newton constant is not positive in general. Indeed, 
this is clear from eq. (12.71) . In fact this was already noticed by Adler in his review [35]. This 
may look surprising that the theory does not automatically predict the correct sign for the 
the strength of the most important interaction. Curiously enough this sign is not fixed by the 
general consistency requirements such as unitarity. Thus, the positivity, if required, may be 
only achieved by an appropriate choice of the fields multiplet. 

2. Even though 1/G(e) is not positive it can be presented as a difference of two apparently 
positive quantities. The Erst one is the Erst term in (12.71) . It represents the contribution of the 
physical degrees of freedom to the renormalization of the Newton cosntant. This contribution is 
obviously positive: each physical degree of freedom, no matter bosonic or fermionic, contributes 
positively to 1 /G(e). The second quantity is what we call the “contact term” c^. It comes 
from the trace of X 1 ' 8 ' 1 in the field operator. We do not have a general proof that cW should be 
positive. However, the simple inspection of the examples considered shows that in these cases 
cis positive. The only possible exception is the scalar field for which the sign of the contact 
term depends on the sign of the non-minimal coupling £. In fact, it appears that the negative 
values of £ for the Higgs boson are phenomenologically preferred [42]. 

3. In certain cases 1/G(e) may vanish. This is so in the case of the conformal scalar in d = 4. 
However, the vanishing of 1/G(e) is not a general consequence of the conformal symmetry. The 
only other conformal field theory we know with same property in four dimensions is the Af = 4 
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SU(N) super-Yang-Mills theory. The respective matter multiplet contains 6 conformal scalars, 
2 Dirac fermions and 1 vector held. 

4. The observable division of all fields in the Universe on the matter fields and those that mediate 
their interactions (gauge bosons, gravitons) is particularly sharp when we discuss their respective 
contributions to the Newton constant. Only the fields-mediators have the non-vanishing contact 
terms and thus produce the negative contributions to the Newton constant. The Higgs boson 
might have some non-trivial non-minimal coupling to gravity and in this respect it might be 
similar to the mediators. 

5. One might worry that the contact term (15.28[) due to graviton is quadratic function of 
dimension d compared to the one due to a gauge vector held which does not depend on d. 
However, the dependence of the contact terms on d should be compared with such a dependence 
of the number of the off-shell degrees of freedom. If we do so we find that both s = 2 and s = 1 
cases follow the same pattern: 

D s (d ) — c^ s \d) = d — s , s — 1,2, (6.1) 

where D s=1 (d ) = d and D s=2 (d) = 'M±ll i s the number of components of symmetric tensor of 
rank 2 . 

6. In the Einstein-Maxwell theory the vector and gravitational perturbations do not decouple 
(this case was considered by Sen in j2Tj, where the main focus was made to the logarithmic 
terms). The relevant second order operator is not diagonal and the result for the effective 
action does not simply factorize on the sum of the contributions due to spins s = 1 and s = 2 . 
Nevertheless we expect that the factorization takes place for the Newton constant and the area 
law in the entropy. The details, however, would be interesting to work out. 

7. An important question is whether the contact terms are gauge independent. The analysis 
given in 5th ref. in [23] shows that in the case of gauge vector fields (s = 1) the contact term does 
not depend on the choice of the gauge fixing term. However, for gravitons, s = 2, it is known 
that the A-term in the effective action may depend on the gauge, see for instance [36], [37] 
(see also [38]). Since the part due to the number of the on-shell degrees of freedom is gauge 
independent this means that the contact term due to spin s = 2 may depend on the gauge. 

7 Entropy 

Let us now discuss the entropy. It is defined for a co-dimension two surface £ by applying 
the so-called replica method. The main case we are interested in is when the surface £ is a 
black hole horizon. Therefore, for simplicity we shall assume that in the curved spacetime there 
exists a timelike Killing vector d T such that it is null at the horizon surface £. The various 
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generalizations are of course possible: when the surface S is a generic surface (not a fixed 
point of any isometry) characterized by some extrinsic curvature, when the Killing symmetry 
is only approximative, valid in a small vicinity of £ and many others. In the replica method 
one identifies the time coordinate r with period 27 tq: and computes the path integral for a held 
in question on a manifold with this new periodicity. The respective quantum effective action is 
function of a, W(a). We are in particularly interested in the case when a — n is an integer. 
The partition function then can be identified with a trace of n-th power of the reduced density 
matrix, 

e -W(a=n) = Trp n (7.1) 

The entropy then is defined by analytically extending (17.11) to non-integer values of a and 
computing the derivative with respect to a, 


S(a ) = ( ad a — 1)IK(o:)| q= i . 


(7.2) 


For integer a = n > 1 the Euclidean manifold in question is an n-fold cover of the original 
spacetime defined at n — 1. It has a conical singularity with angle deficit 5 = 27t(1 — a) at the 
surface £. Only this singularity contributes to the entropy (17. 2 |) . 

Technically it is more convenient to work with the case when a = 1/N < 1. The manifold in 
this case is an orbifold. Namely this space arises when we consider the gravitational held created 
by a cosmic string of positive tension. Geometrically and topologically these two manifolds, the 
orbifold and n-fold cover, are quite different. The conical singularities however in these two 
cases are analytically related and to some extent one may not worry about the differences 
between them. The attention however is required when the higher spins are considered. Then 
some non-analyticity in passing from a < 1 to a > 1 may arise. We shall discuss this issue 
later in the paper. 

As we have seen the partition function for a held of spin s reduces to calculation of deter¬ 
minant of some Laplace type operators 


-\/ 2 s ab + x^ b , 

(s) 

where the potential term X AB is constructed from the Rieamann tensor, 


\d s ) _ D 

^AB a 


AB 


(7.3) 


(7.4) 


Thus, the calculation of the effective action (17.11) and the entropy (17.21) reduces to calculation of 
the heat kernel for the operator of the type (j7.3[) on a manifold with a conical singularity. The 
main question here is how to treat the potential term (17.4[) . On a conical manifold the Riemann 
curvature has a delta-like contribution due to the singularity m 


= 2tt(1 - a)((n a ,n^)(np,n u ) - (n Q ,n I/ )(n / 9 n / 1 ))5 E + R r ^„ , (7.5) 
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where (n a , np) = n“, a = 1,2 is a pair of vectors normal to E, is the regular part 

of the curvature which remains when a = 1. 

The question now is whether we should take into account this singular term in the curvature 
(17.5p when compute the effective action and the entropy or this term should be ignored? There 
are two possible answers to this question and, respectively, two different entropies can be defined. 

7.1 Entanglement entropy 

One possibility is to treat the potential term (17.4)1 as any other potential term. This suggests 
the following prescription. First we consider operator (17.3ft with some unspecified potential 

D = -V 2 5 AB + V AB , (7.6) 

compute the partition function Z(a) for this operator on a conical space, deduce the entropy 
S = — ( ad a — 1) In Z(a)\ a= i and only at the end make the substitution: V AB = X AB . In this 
way we completely ignore the singular term in (17. 4p due to the singularity in the Riemann 
tensor (17.51) . This prescription is motivated by that simple observation that the construction 
of the reduced density matrix from a vacuum state in Minkowski space-time does not appear 
to be affected by the presence of the Riemann tensor in the potential term (17.4ft . Indeed in 
flat spacetime this potential terms is vanishing and the reduced density matrix appears to be 
exactly the same as for an operator with X^ s > = 0. The entropy defined by this prescription we 
shall call the “entanglement entropy”. It corresponds to a well-defined density matrix and, in 
particular, this entropy is positive definite as we show below. 

The calculation of the partition function Z(a) and the entropy S(a) can be performed using 
the heat kernel technique. The coefficients in the expansion (12.31) . being considered on a conical 
space, are modified by certain surface terms, 

a n (a) = aa n (a = 1) + a^(a), (7.7) 

where a n (a = 1) are the coefficients in the heat kernel expansion on regular spacetime and 
a^(a) are the surface terms given by integrals over the entangling surface E. In this paper, 
for the analysis of the leading divergences of the entropy, we shall need only the first n — 1 
coefficient. For an operator of the type (17.6[) we have that 

of (a) = [ 1, (7.8) 

3cc is 

where D(d ) = S AA is the dimension of the representation. If the operator (17.6)1 is acting on 
fields of half-integer spin then the surface term (17.8[) is negative compared to the case of fields 
of integer spin, this behavior is taken care by the factor (—1) F in (17.8ft . Notice that the first 
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coefficient (17.81) does not depend on the potential V. The case of the Dirac fermions was first 
considered in ESI- Here we closely follow the discussion in |39] . 

It should be noted that for spins s = 2 and s = 3/2 there have been found in [39] (see 
also j40j ) some extra surface terms additional to (17.8[) which do not disappear when a = 1. 
These non-analytic terms are not written in (17.81) . They will be discussed later in section 17131 
where we argue why these terms should not be present in the n-fold cover construction, i.e. 
when a — n > 1. 

In the examples considered above the complete partition function for a held of spin s contains 
a product of positive and negative powers of determinants of the Laplace type operators D*, 

(CD, 

Z (s) = IJ det^ (fc)/2 D fe , (7.9) 

k 

where u(k) are some numbers, positive or negative. For each operator one has the heat 
kernel expansion (17.71) . (17.81) . The leading UV divergence of the surface part of the effective 
action then can be computed as 

L 1 • (7 ' 10) 

where Dk(d) is the dimension of the representation in the operator D^.. If the operator is 

acting on fields of half-integer spin its contribution in (17.101) is counted with negative sign. In 
all cases considered in this paper we have observed that the sum 


£>*(-! ) F D k (d) = V s (d) (7.11) 


is the number of on-shell degrees of freedom of spin s. Using this property and the form (17.101) 
for the surface action we find the leading term in the entanglement entropy due to spin s, 

1 1 V 8 {d) 


c( s ) _ 

°ent 


(47r)V(d-2) 


2 


6 


-H(E), 


(7.12) 


H(E) is the area of surface £. We see that only physical degrees of freedom contribute to this 
entropy. Moreover, the entropy (17.12)1 is obviously positive as it should be provided the entropy 
is determined by a well-defined density matrix. The other observation is that the entanglement 
entropy (17.121) is gauge independent. 

The generalization to the massive fields is straightforward and is carried by the multiplication 

2 

of the heat kernel by a factor e~ m T . The complete result for the entanglement entropy due to 
a set of spin s and mass m s is 


Sent — 


1 'E > s,m(d) 

2 ^ 6 



dr 


(47rr) d / 2 


44(E). 


(7.13) 
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Comparing the area law (17.1211 . (17.1311 and the renormalization of the Newton constant (12.711 . 
(13.191) we see that if there are no fields in the multiplet for which the contact terms are 
non-vanishing then the entropy (j7. 12[) takes the Bekenstein-Hawking form 


S 


ent 


MX) 
4 G(e) 


(7.14) 


expressed in terms of the induced Newton constant G(e). This is in fact a part of the theorem 
proved in P2] about the simultaneous renormalization of entanglement entropy and the Newton 
constant. In the presence of the gauge fields (s = 1 and s = 2) the situation, however, is 
different: due to the contact terms one can not in general write entanglement entropy (17.12)1 in 
the form (17.14ft . 


7.2 Conical entropy 

The other possibility is to take into account the singularity in the curvature (17.5jl even in flat 
spacetime. This reduces the problem to considering the quantum held theory coupled to a 
delta-like potential. The complete treatment of the problem might be quite complicated (see 
for instance HP). We are however interested in the case when the angle deficit is very small 
so that the singular term in (17.51) can be treated as a perturbation. For a non-minimal scalar 
this approach was suggested in [T9]. The perturbation theory for a delta-like potential is not 
well defined since the higher powers of the “perturbation” are singular non-iutegrable functions. 
This problem however is absent for the terms linear in the perturbation so that to linear order 
in (1 — a) the perturbation theory gives a meaningful result. The entropy which appears in this 
calculation we shall call the conical entropy since it is defined as a response to introduction of 
a small angle deficit in the effective action. The conical entropy is not related to a well-defined 
reduced density matrix and it does not have, in general, a clear statistical meaning. It has 
however other advantages which will be discussed below. 

It should be noted that for the fermionic fields (s = 1/2, 3/2) the second order operator 
is obtained by squaring a first order Dirac type operator for which there is no a non-minimal 
coupling to gravity. Therefore, it appears reasonable that the singular terms in the curvature 
should not be taken into account in the corresponding second order operators. Thus, below in 
this section we shall consider only bosonic operators. 

Considering the singular term in (1 7.5 j) as a perturbation to linear order we find a modification 
of the surface terms in the heat kernel coefficients of the operator (17.31) . (17.41) . This modification 
arises in exactly same way as in the case of the scalar field considered in [19]. Generalizing this 
consideration to a more general potential term (17.41) and taking into account (17.51) we find that 
to leading order in (1 — a ) 

a i = a l( 0 ) - 47r ( 1 - «) [ Tr U a ^(n a , n^) (rty, n„) , (7.15) 

Jt, 
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where af( 0) is the unperturbed coefficient fl7.8D . and the trace is taken with respect to indexes 
AB of matrix U. Now we notice that the trace of the potential term (17.41) reduces to a Ricci 
scalar, 


TrX (s) = TiU a ^ u R aPflp = cR. 

For the Ricci scalar the formula (17.5ft reduces to 

R = 4vr(l - a)S s + R reg . 

Applying (17.5p . (I7.17|) to both sides of (17.16[) we fold that 

f Tr U aPllu (n a ,n li )(np,n l/ ) = c [ 1 
J s Jr, 

and the surface term in the heat kernel coefficient for a bosonic operator (17.31) reads 


r _ /V(l-a 2 ) 

Lt-I - l 


i - | 77--- D,(d) - 4jt(1 - a)cj j 1. 


(7.16) 


(7.17) 


(7.18) 


(7.19) 


Let us consider some examples. 

1. Scalar operator —V 2 + C.R. This case was already considered in [TU] . In this case c = £ 
and we have for the surface coefficient 

«? = J ^ ~ 4tt(1 - ot)£ \ I 1 = 4tt(1 - a)(^m f 1 + 0(1 - a) 2 . (7.20) 


2. Vector Laplace operator Ai = —'V 2 5 fll/ + R^ u . In this case c = 1 and we have that 




(7.21) 


3. Laplace operator on p-forms A p = (d5 + 5d) (14.51) . In this case c = C^__\ and hence 

V (1 — a 2 ) 


a? = 


a 


-cj - 477(1 - ^ 1 


(7.22) 


4. Lichnerowicz operator (15. 3 p acting on symmetric tensors of rank 2: c = (d + 2) (as can 
be seen from (15.261) ) and 

V (1 — a 2 ) d(d + 1) 


O] = 


a 


— 47t(1 — a)(d + 2) ) j 1. 


(7.23) 


These results for the vector and Lichnerowicz operators should be compared with the earlier 
results obtained by Fursaev and Micle [39] (and re-derived recently using some other methods 
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in [IQ]). The analysis used by these authors is based on exact spectrum of the operators and 
thus is supposed to be exact while our analysis above is based on a perturbative treatment of 
the singularity. The comparison shows that these results for the analytic terms (vanishing when 
a = 1) completely agree. This agreement for the Lichnerowicz operator however is not extended 
to the non-analytic term (independent of a and non-vanishing when a = 1) found in [[39]. We 
discuss the non-analytic terms in the next section. 

We can now compute the surface terms in the effective action of the fields of spin s considered 
earlier in the paper. The replica formula (17.2p applied to the surface effective action that 
corresponds to the surface heat kernel coefficients gives the following result for the leading 
area term due to a field of spin s (where we have to take into account all contributions due 
non-minimal coupling in the ghosts operators): 



1 1 

(47r)^(<i — 2) ed ~ 2 



-c (s \d)^j 71(E). 


(7.24) 


Including the fields of various spin s and mass m s we obtain 


q 

mn. 



V S ,m{d) 

6 



dr 


(47j- r )d/2 


-4(E) 


(7.25) 


for the complete conical entropy. We see that this expression repeats exactly the renormalization 
of the Newton constant (13.191) so that we can write it in the form 


MX) 

■ 4 G(e) ' 


(7.26) 


Namely this property of the conical entropy makes it so attractive. We, however, notice immedi¬ 
ately that the conical entropy (17.241) is not positive definite and thus likely can not be obtained 
from a well-defined density matrix. The other observation is that the conical entropy (I7.24p . 
(I7.25P is gauge dependent since the contact term for graviton may depend on the gauge. 


7.3 Puzzle of non-analytic terms: orbifolds vs n-fold covers 

As we have already mentioned, Fursaev and Miele in their 1996 paper have found some non- 
analytic terms in the heat kernel of operator for spin 3/2 and of the Lichnerowicz operator [39]. 
The similar results were recently obtained in [40j using a somewhat simpler method that can 
be further extended to compute the entropy in string theory. 

The non-analytic term in the heat kernel for fields of spin 3/2 is proportional to a and thus 
does not produce any contribution to the entropy. Therefore, we shall discuss only the case of 
the Lichnerowicz operator for a field of spin s = 2. 

Let us briefly outline the way this term emerges in the analysis of [39]. They consider 
the two dimensional case as the basic one and use the heat kernel obtained in two dimensions 


22 

















to derive the coefficients in higher dimensions by considering the products of two-dimensional 
space and (d — 2)-dimensional surface £. This approach is sufficient if we are interested only 
in coefficient a\. Thus, they consider a singular two-dimensional sphere which is obtained 
from the regular sphere with coordinates (9,t) by identifying angle r with period 27tcx. They 
specifically consider the case when cx < 1 that corresponds to an orbifold. On such a singular 
sphere Fursaev and Miele construct two scalar functions 

= sin" 9e ±iaT (7.27) 


that are eigen-functions of the Laplace operator 


V 2 <£ = 


(1 + cx) 


cr 




(7.28) 


From these two functions one may construct 4 vectors ( V±)n = e^d v (j>^ and (V^) M = ■ 

For regular sphere (a = 1) 2 vectors V ± are the Killing vectors and the other two vectors V ± 
are the conformal Killing vectors. Therefore, they can not be used to produce any physical 
gravitational perturbation. On the other hand, when a / 1 these 4 vectors are no more the 
Killing vectors and thus can be used to construct 4 gravitational modes 


(ftJU = v„(r±)„ + v„(v±)„, (hi v 


V„(V±)„ + ?)“ • 


(7.29) 


These modes are new compared to those gravitational modes which already existed for regular 
sphere (cx = 1). Each mode contributes 27 t to the heat kernel coefficient cq. For a higher 
dimensional space which is product of two dimensional space and a (d — 2)-dimensional surface 
£ this results to modification 

af = + 87t J 1 , (7.30) 

where af L is the surface coefficient (17.231) . We see that this new term does not depend on cx 
and does not disappear when cx = 1 . 

It is an important question whether these new modes (17.29[) are normalizable. One finds for 
these modes that hee ~ (sin 0)^^ 4 so that they are normalizable with measure sin 6d6 only if 
cx < 1. The latter condition is satisfied when S\ is orbifold. On the other hand, these modes 
are non-normalizable if a > 1 and in particular they should be excluded when ex = n is an 
integer. Thus, for an ?r-fold cover of the two-dimensional sphere the gravitational modes are 
the same as for the regular sphere. The authors of [39j have speculated that for a > 1 there 
may exist a self-adjoint extension of the Lichnerowicz operator for which the heat kernel is 
analytically continued to all values of a. However, it is more reasonable to think that, once 
we are dealing with some non-analyticity (between cx < 1 and cx = 0 cases) anyway, the case 
of cx > 1 is analytically related to the regular case cx = 1. This is achieved by excluding the 
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modes (j7.29j) from the spectrum that results in the same heat kernel (j7.23jl as we had before, i.e. 
without the non-analytic term 87 T. This consideration demonstrates that the orbifold and an 
n-fold cover are not in general analytically related to each other. In the context of the entropy 
calculation that uses the replica trick we are dealing with an n-fold cover. We conclude that for 
the entropy calculation the non-analytic terms should not appear. The absence of these terms 
is actually important for the renormalization of the entropy, both entanglement and conical. If 
the non-analytical terms were present neither of these entropies would be renormalized in the 
same way as the Newton constant. 


8 Can the Bekenstein-Hawking entropy be equal to en¬ 
tanglement entropy? 


Let us first consider the case when the only fields present are the matter holds, scalars or 
fermions, that couple minimally to gravity. For simplicity we shall consider the massless fields. 
The generalization to the massive held is straightforward. Then, as we have already noted 
earlier, the entanglement entropy and the Newton constant are renormalized in exactly same 
way so that, provided there is no bare Newton constant, i.e. the entire gravitational coupling 
G ren is induced, we have that 

1 _ 1 _ 1 AT V s( d ) /Q 1 \ 

4G ren ( 4 7 t) ^(d — 2) e<i ~ 2 S 6 ’ 1 ; 


where N s is the number of helds of spin s. Here e is the UV cut-off. By relation (18.TJ) it is 
hxed to be proportional to the Planck length. If there is a black hole with horizon E in the 
induced gravitational theory then it has the Bekenstein-Hawking entropy defined in terms of 
the gravitational coupling constant G ren (18.111 


Sbh — 


/1(E) 


( 8 . 2 ) 


where /4(E) is the area of E. On the other hand, the entanglement entropy of E due to the 
same set of helds is exactly the Bekenstein-Hawking entropy 


Sent — 


21(E) 

AG r en 


= S. 


BH 


(8.3) 


This is guaranteed by the fact the Newton constant and the entanglement entropy are renor¬ 
malized in exactly same way. Thus, if there were only matters helds, scalars and fermions, in 
the Universe then the above mechanism would be the natural explanation for the statistical 
origin of the black hole entropy. 

However, in the observable Universe there exist not only matter helds but also the helds which 
mediate the interactions, gauge vector bosons and gravitons. For these helds the renormalization 
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of the Newton constant and the entanglement entropy is different due to the presence of the 
contact terms c^ s \d) that contribute negatively to the Newton constant and do not appear at 
all in the entanglement entropy. The presence of these fields in our Universe makes impossible 
this simple identification of Bekenstein-Hawking entropy and the entanglement entropy that 
was just outlined. At this point it is tempting to suggest (as for instance in P) that instead 
of the entanglement entropy we should consider the conical entropy which, as we discussed this 
above, repeats identically the renormalization of the Newton constant. Therefore, extending 
the identification of the observable Newton constant (18.11) to other spins, s = 1 and s — 2, we 
would still have the equality between the Bekenstein-Hawking entropy and the conical entropy, 

Sbh = S con . (8.4) 


The only problem with this idea, as we already mentioned this earlier, is that the conical entropy 
is not a priori positive definite (the positivity would impose extra constraint on the multiplet 
of fields) so that, most importantly, it is not determined by a well-defined density matrix and 
hence the desired statistical meaning of the entropy is lacking. 

In general, when the bare Newton constant G is present, the renormalized gravitational 
coupling is defined by the relation 


4 G r 


1 

4 G 


(47r)V(d-2) 


cd-2 




T> s (d) 


-c« 


(d)) 


(8.5) 


Notice that in the contact terms we include also the contribution of scalar particles, c^ =0 ^ = £, 
due to the phenomenological possibility [T2] that the Higgs boson may have a negative (and 
even large |£| 3> 1) non-minimal coupling to gravity. 

The question we want now to ask is whether it is possible to modify the above mechanism 
and have at the same time two things: 

1. ) Gravitational constant is defined according to (j8.5[) and, respectively, the Bekenstein- 
Hawking entropy is defined as Sbh = A(E)/4G> en in terms of this constant. 

2. ) Entanglement entropy defined as 


Sent — 


(4vr)^(d-2) 


c d -2 




( 8 . 6 ) 


is identical to the Bekenstein-Hawking entropy, 


Sbh = S, 


ent 


(8.7) 


The answer to this question is affirmative provided the following consistency condition 


1 

4 G 


1 1 

( 47 t )^ {d — 2 ) e<1 ~ 2 


Y, N ,c M (d) 


( 8 . 8 ) 
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is satisfied. Since the bare Newton constant G is a free parameter in the theory it always can 
be chosen to satisfy the condition (18.8b This condition simply means that the bare Newton 
constant G cancels completely the contribution due to the contact terms in (18.5p so that the 
renormalized Newton constant is defined only by the physical degrees of freedom 


1 = _ 1 _\ ^ 'Dsi.d) 

4 G ren (47r)^(d — 2) e<l ~ 2 S 6 


(8.9) 


in the same way as in ( 18.lj) . The ratio of the bare and the renormalized constants is determined 
only by the field multiplet 


Gren _ Es^Kd) 

G 52 s N s V s (d) ' 


( 8 . 10 ) 


We see that the presence of the bare Newton constant G is important for this mechanism to 
work. In the complete Bekenstein-Hawking entropy there are two pieces which do not have a 
statistical meaning: the one due to the bare constant G and the other due to the contact terms. 
With the condition (18.8p these two pieces mutually cancel and what is left has a perfect statistical 
meaning as entropy of entanglement. As we have already mentioned above, the generalization 
of this mechanism and, in particular, of relations fl8.8p ~ (18. 101) to the case of massive fields is 
straightforward and does not bring any new principle issues. 

A nice bonus of the proposed mechanism is that both the entropy (18.61) and the observable 
Newton constant (18.91) are determined only by the physical degrees of freedom and, thus, do 
not depend on the gauge. 


9 Conclusions 

As we show in this paper the analysis of the Newton constant and the entanglement entropy 
can be done in terms of the fields which appear in the low energy description. The whole 
picture thus does not seem to depend on the details of the fundamental theory. In fact, this was 
anticipated by many authors who noticed the universality of the Bekenstein-Hawking entropy 
which does not seem to depend on how the low energy gravitational physics is completed in the 
higher energy regime. In the present context, however, the fundamental theory should give us 
a concrete realization of the mechanism proposed in this paper. In particular, we expect that 
the low energy parameters such as the bare Newton constant G will arise unambiguously in the 
fundamental theory and hence the relation (18.81) should be reproduced in a concrete calculation. 

The other lesson which we should learn from the above consideration is that the orbifold and 
the rt-fold cover constructions are not always analytically related. Even though the orbifolds 
look easier to make calculations on a better understanding of the n-fold covers is probably 
needed to make a further progress. 
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